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Abstract 

We study the pair creation of point -particles and strings in a time-dependent, weak gravitational 
field. We find that, for massive string states, there are surprising and significant differences between 
the string and point-particle results. Central to our approach is the fact that a weakly curved 
spacetime can be represented by a coherent state of gravitons, and therefore we employ standard 
techniques in string perturbation theory. String and point -particle pairs are created through tree- 
level interactions between the background gravitons. In particular, we focus on the production of 
excited string states and perform explicit calculations of the production of a set of string states 
of arbitrary excitation level. The differences between the string and point-particle results may 
contain important lessons for the pair production of strings in the strong gravitational fields of 
interest in cosmology and black hole physics. 
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I. INTRODUCTION 



The aim of this paper is a calculation the pair production of strings in a time-dependent 
background. This phenomenon is of paramount importance in cosmology, particularly during 
(p)reheating after inflation, and near cosmological singularities. Usually one studies pair 
production using an effective field theory approach. Here we study a spacetime with weak 
gravitational fields and compute the pair production of excited string states using worldsheet 
methods. This "stringy" calculation gives significantly different results than those found in 
the effective field theory approximation. Our results suggest specific signatures of the string 
pair creationphenomenon that may survive into the regime of strong gravitational fields. 

Recently [l|, AS the pair creation of strings in a cosmological setting was estimated using 
effective field theory methods. In this work the excited states of the string are treated as 
massive particles for which the usual field theory calculation applies. While the production 
rate of a single massive state is small, the exponential Hagedorn growth of the number of 
string states can give rise to a significant cumulative effect. Here we test this effective field 
theory approximation in a weak-field limit, where string pair creation may be computed 
exactly. Although our real interest is strong gravitational fields, such as those that arise in 
cosmology, the formalism does not yet exist to calculate pair production of strings here (see 
^| for an interesting approach to this and for a calculation in a more specific background). 
Nevertheless interesting surprises can be seen already for weak fields. Among them, we find 
that the enhancement of the production of excited string states due to the growing density 
of states is greatly suppressed due to the very mild hard scattering behavior of the string, 
which in turn arises from the fact that the string is an extended object. 

To carry out our calculation of string pair production, we use the familiar idea that 
a curved geometry can be represented coherent^ of graves g J FOT we ak 
gravitational fields this gives a prescription relating S'-matrix amplitudes on a weakly curved 
spacetime to a sum of those on Minkowski spacetime with multiple graviton vertex operators 
inserted. To lowest order in the closed string coupling, the pair production of strings is thus 
the four-point process gg — > AB, with two gravitons g interacting to produce two strings 
in states A and B. This approach has the advantage of leading to a controlled calculation 
using standard string perturbation theory, but there are many other conceptual and technical 
obstacles to an understanding of strings in general time-dependent backgrounds. There are 
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subtleties involved in interpreting the S'-matrix for strings on curved backgrounds (see, for 
example, jsl 0]). More generally, in any quantum theory including gravity the inevitable 
formation of singularities as predicted by the singularity theorems [llOlJlll] implies the absence 
of asymptotically trivial "in" and "out" regions required to define an S-matrix. We construct 
a specific geometry describing the collision of two small amplitude plane waves which we 
embed in an exact nonlinear solution of Einstein's equations. An analysis of the full solution 
shows the formation of a spacelike singularity in the future of the collision of the plane 
waves, and also how by making the amplitude of both waves sufficiently small we may push 
the singularity arbitrarily far into the future of the initial collision region. Consequently 
perturbation theory is valid for an arbitrarily long period, and we believe the S'-matrix 
formalism can be used to understand the physics in this region. This is as it should be, for 
the possible formation of a big crunch in our future should not prevent us from understanding 
physics today! 

In this work we primarily focus on the pair production of long-wavelength string states. 



Since the theory of pair creation in field theory is well developed [12|, Q], and since string 
theory should have an effective field theory description valid for A ^> Va', one might ask 
why we expect there to be differences between string and point particle pair creation in this 
regime. When viewed as the four-point process gg — > AB, then certainly for massless states 
A and B, with wavelengths A ^> \fa! ', all momentum scales in the scattering process are well 
below the string scale. On the other hand, to produce excited string states, with masses 
m ~ then the typical graviton momenta and momentum exchange will in fact be 

greater than 1 / \fa' ; . This is beyond the validity of effective field theory, and indeed previous 
investigations have revealed uniquely "stringy" behavior in this regime 15, Thus, 
it is reasonable to expect the pair production of massive strings might be fundamentally 
different from the pair production of massive particles. 

Our calculation reveals two specific differences between the string and point-particle 
results. As we have mentioned above, and will discuss in more detail in Section I1V| the 
pair production of strings is suppressed due to the mild hard scattering behavior of string 
amplitudes. Another intriguing difference relates to the production of identical string states. 
At tree level in field theory only identical particle pairs are produced. A surprising feature 
of the string case is that production of pairs of strings at the same excitation level are 
suppressed, going to zero as A — > oo. Instead, most created pairs are of different mass levels. 
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Both of these features can be explained in the context of certain simple features of the string 
spectrum and scattering amplitudes, as we discuss in more detail below. 

This paper is organised as follows: In Section |H] we discuss the perturbative approach 
to calculating pair production, using field theory as a guide to the string calculation. In 
Section 11111 we discuss a specific example of a background, describing the collision of two 
plane waves, for which the perturbative calculation may be applied. 

In Section llVl we perform the explicit string calculation and discuss its properties. Finally, 
we present our conclusions in Section El We also include two appendices with additional 
details on some elements of our work. The four-point amplitude describing pair produc- 
tion is discussed in Appendix |XJ and a more detailed account of perturbatively computing 
Bogolubov coefficients is given in Appendix IH1 



II. PAIR PRODUCTION IN WEAK GRAVITATIONAL FIELDS 



In this section we discuss how pair production may be studied perturbatively. Certainly 
there are examples where the production rate is nonperturbative; these include Schwinger's 
classic calculation of e + e~ production in an electric field and more recent calculations 



with open strings [18, 
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20j . An example with a different behavior is provided by the pair 
production of scalars of mass m in a flat FRW spacetime that undergoes a brief period of 
expansion j^jj. In this example we have the metric 



ds 2 = a(i]) 2 {—drj 2 + dx 3 ) , a 2 (r]) = I + bt&nh(pr]), 

and the number of particles produced in each mode k = (u, k) is found to be 

sinh 2 [tt(uj +00 - ^_ oc )/2p] 



(1) 



Nt 



sinh(7rcj +00 /p) sinh(7ru;_ 00 /p) 



where 



u}(rj) = \k\ 2 + a(r]) 2 m 2 



1 1/2 



U±c 



uj(±oo). 



(2) 



(3) 



For b <C 1 the number of pairs produced goes like ~ b 2 m 4 for small b. Thus there are 
backgrounds in which is is plausable that pair creation could be studied within a perturbative 
framework. 
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A. Field Theory 



In the case of field theory, studying pair production can be realized quite straightforwardly 
through the addition of new vertices in the Feynman diagrams of the theory. This point 
of view is developed in more detail in Appendix [HI where it is shown that this technique 
yields precisely the same results as the more standard Bogolubov coefficient calculation. 
The field theory case is useful for our present purposes as it provides a prototype for the 
string calculation. One expands the metric about Minkowski space 

g^u = + h$ + h$ + • • • , (4) 

where h^l is the first order metric perturbation, and nfi} + . . . , the corrections required at 
successively higher order so that the full metric g^ satisfies the Einstein equations. If we 
consider a massive minimally coupled scalar field with action 

S ~P = -\j (^<R0 + m 2 2 ) y/=g d D x, (5) 

Then the expansion (£Q) of the metric gives a corresponding perturbative expansion of the 
action, 

S, = Sf + + ••• (6) 

as well as a series of terms at n th order in the perturbation. 

This allows a perturbative evaluation of the n-point functions in the background, de- 
noted by (■ • • )b, in terms of the same quantities (■ ■ • )o in Minkowski space with additional 
insertions. For any combination of fields O, the expansion of the action © inserted into 
the path integral yields, 

(0) b = (O) + J {TO k ^%^ + j (TO„Y®) ^ + ■■■ (7) 

where and E^ 2 ^ are the new vertices that appear to first and second order in perturba- 
tions. These are determined by the expansion (JH| of the action 

£«=i<S«, (8a) 

(8b) 

The statements that we have made for the n-point functions have analogues for the S- 
matrix. However, defining an S-matrix involves a choice of asymptotic "in" and "out" 
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£( 2 > = iS 



(2) 



s 



(1) 



states. Therefore this prescription really only makes sense when the spacetime in question 
is asymptotically Minkowski. In other cases it may be necessary to modify the definition of 
asymptotic states. 

At this point, it is possible to see how even a weak gravitational background can lead 
to pair production. Pair production is the process with no "in" particles and two "out" 
particles. While forbidden by momentum conservation in Minkowski space, in the perturbed 
spacetime each term in the ^-matrix carries three momenta; the two momenta of the "out" 
particles and the momentum carried by the new vertices E^. Thus two positive energy 
particles can appear in the "out" state if one of the terms Ej. supplies the "missing" 
momentum. So, the amplitude to pair produce two particles with momenta ki,k 2 will be 
nonzero provided that we have 7^ 0. More physically the missing momentum is 

coming from the background gravitons. 



B. String Theory 

Following the same procedure as the field theory calculation, one may insert the pertur- 
bative expansion of the metric g^ v = rj^ + hfy + h$ H into the conformal gauge string 



worldsheet action 



S = -— I gixv {X)dX»dX v d 2 z. (9) 



Now, the successive metric perturbations hffl are chosen to satisfy the string beta function 
equations. 

As is well known ((J Q] the expansion of the metric leads to the insertion of graviton 
vertex operators in worldsheet correlation functions. For example, if we wish to study the 
S'-matrix element involving two string states A and B then we must compute the worldsheet 
correlation function 

1 \ d D k' 



o (2 



V A (k A )V B (k B )^ b = (v A (k A )V B (k B ))+^- J h$(k')(v A (k A )V B (k B )V^(k'] 

^ (10) 
where V A)B {k AtB ) are the vertex operators for the created strings in states A and B, V£ u (k') 
is a graviton vertex operator integrated over the worldsheet representing the metric pertur- 
bation, and we have absorbed a factor of g c in each of the vertex operators. The metric 
perturbation h$ appears as the polarization tensor of the background graviton. 
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The first-order contribution to pair creation (|10|) vanishes when A and B are massive, due 
to momentum conservation when the graviton is on-shell (as it must be for a consistent string 
amplitude). Nevertheless, (|1U|) is related to the phenomenon of "particle transmutation," by 
which a string can change its mass and spin as it passes through a gravitational background 
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23, 
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25| . We will not investigate this phenomenon in the current work, but for weakly 
curved spacetimes it could in principle be studied within the S"-matrix framework we apply 
herein. 

Ideally we would be able to continue this expansion to higher order in the metric per- 
turbation, by analogy to the field theory case. However, at second order issues arise in the 
string theory that make the situation somewhat more complex. For one, consistency of the 
of the string amplitude (jl(Jj) requires that each of the vertex operators appearing in be a 
conformal tensor of weight (h, h) = (1, 1), or in other words satisfy the physical state condi- 
tions. As is well known, this is equivalent to imposing the linearized Einstein equations on 
h$- At second order in the expansion of worldsheet amplitudes, a graviton vertex operator 

(2) (2) 

with polarization h)J appears. By construction, the second-order perturbation h^l satisfies 

(2) 

the second-order RG equations. In general, this means that h)il with not in fact satisfy 
the physical state conditions. Thus, the string amplitude with the corresponding vertex 
operator will not be consistent. 

Another issue that arises is the presence of a nonzero "tadpole" for some fields. Taking 
our prescription literally, to compute the S'-matrix elements on our background to second 
order involves inserting two graviton vertex operators into the relevant correlation functions. 
For massive fields A, the tadpole is zero at first order due to momentum conservation. At 
second order in the metric perturbation, the tadpole for some state A in the background is 
therefore given by 

(v A (k A )) b = (V j4 (A;a)V 9 (A; 1 )V 9 (A; 2 )) + • • • (11) 

with V g (k) = ^—h^uV^ u (k). This expression does not vanish for the examples we consider. 
A nonvanishing tadpole is an indication that the classical equations of motion are not exactly 
obeyed by the background solution. This is to be expected, as here we will use a background 
solution that solves the leading-order equations of motion for massless string excitations, 
but we have ignored the massive excitations and higher corrections. The presence of the 
tadpole is telling us that the equations of motion for these massive fields forbid setting them 
to zero throughout spacetime. Furthermore, in a curved background the vertex operators 



and states should change, which we have not taken into account. We expect that the failure 
of the tadpole to vanish is due to some combination of all of these effects. 

Given these new complexities, it appears that problems are beginning to develop with the 
technique of inserting new vertex operators at second order. Nevertheless, by considering 
each of these new problems, we will argue in Sections 111 CI that useful results can be obtained 
from the second order expansion. Therefore, for the purpose of studying pair production we 
will use the prescription 



where V^bC^s) are the vertex operators for the pair created string states A and B. In 
defining this S-matrix element it is understood as usual that the positions of three of the 
vertex operators are fixed and the remaining one is integrated over the worldsheet. This 
prescription does not involve a term at first order in the metric perturbation, since the 
pair production of massive states is forbidden to this order by momentum conservation. 
Furthermore corrections to the vertex operators and background solution, which should be 
included in a more complete analysis, only affect this S-matrix element at higher orders in 
h {1) 

C. Coherent states 

The existence of non-zero "tadpoles," or 3-point functions (V s V 9 Vat) where Vn denote the 
vertex operators for level N string states, implies that turning on a nontrivial gravitational 
background will, at second order, produce a nontrivial background of excited states. The 
question we must address is whether we can still trust the pair production results that we 
will extract at this order. By studying an analogous situation in field theory, we will argue 
that the corrections to the pair production amplitudes appear at higher order, and thus our 
second order results may be trusted. 

The production of excited string states can be modeled in field theory as the pair produc- 
tion of a massive scalar field x Yia - interactions with a massless scalar field H, representing 
the graviton. This can be realized with the action 



(AB\S\vac) 



(47T£ C ) 



1 




d D k' d D k" 
o(27r) D (27r) D 



(12) 



S 




l -m\ 2 + XH\ - \^H\ 2 



) 



(13) 



S 



The presence of the A coupling implies that we cannot turn on a background field for H 
without also exciting a background value for x, as a consequence of the equations of motion 



-DH = 2\H X +■■■ (14) 
-n X + m 2 X = XH 2 + (15) 

The \i coupling gives rise to pair creation of x quanta in the presence of a nontrivial 'gravi- 
ton' background H. Before calculating this effect we must first correctly deal with the 
background. Assuming x — + 0(H 2 ) these equations may be solved perturbatively to 
second order in H as (here J dk = J d d ~ 1 k/(2uj), uj = \k°\ which is valid for massless and 
massive particles) 

H(x) = H (x) = J dk (a*(k)e ik - x + a(k)e~ ik - x ) + 0(H$), (16) 
X (x) = \f d d yG ret (x,y)H 2 (y) + O(H 3 ), (17) 

where G ret (x,y) is the retarded propagator for a scalar of mass m and we have chosen the 
following boundary conditions 

lim x(x) = 0. (18) 

It is well known that in field theory a nontrivial background may be described by a coherent 
state. The initial state contains only gravitons and may be described by 

\vac,t^ -oo> = e -U*My* Qot |0). (19) 

Whereas the final state is a coherent state of both 'gravitons' and x quanta 

\vac,t^ +00) = e -hfdk\P\ 2 e -IdH\a\ 2 e fdkaa^ e fdkl3d^ y ^ 

Here a^/a, d^/d are the conventionally defined creation/annihilation operators for H and x 
respectively, and (3{h) is the generated background for x 

lim x (x) = [ dk (p*(k)e lk - x + f3{k)e~ ik - x ) . (21) 

t^+oo J 

We now compute the amplitude to create a pair of x quanta. Because of the nontrivial 
background for X as t — > +oo we must redefine the creation/annihilation operators as 
d = d — p. This is the analogue of cancelling the tadpole in the string theory. The (3 are 
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also determined by the condition that the amplitude to create a single particle must vanish 
(tadpole cancellation) 

(vac,t — > + oo | V2ujdS \ vac, t — ► — oo) = 0. (22) 

Now the amplitude to create a single pair of \ quanta with momenta (k 3 , k^) is 

A[t>ac — > 2] = (vac,t — > +oo| V2u;3d(/c3)v / 2w 4 d(/c4)S'|i>ac, £ — > — oo), (23) 

where S is the S'-matrix. Utilizing the fact that the amplitude to create a single particle must 
vanish and expanding A[vac — *> 2] to second order we find after making uses of conservation 
of energy and momentum 

A[vac^2]= J dk\ J dk 2 a(fc 1 )a(fc 2 )(0| v / 2^rf(A;3)v / 2^rf(A;4)^ t (A;i)a t (A; 2 )|0) + 0(H$). 

(24) 

This is a very simple result: to second order in perturbations, the amplitude to create 
a pair of particles is determined by the amplitude for two gravitons to convert to two x 
quanta. Furthermore, at second order the result in independent of the modifications to the 
creation operators required to cancel the tadpole. This suggests that our string results will 
be reliable, even though we will not carry out the modifications to the vertex operators 
required to cancel the tadpole in the string case. 



III. BACKGROUND 



In this section we shall construct a specific example of a perturbative geometry for which 
our approach may be applied. In doing so we shall also uncover some of the inevitable 
limitations of the perturbative approach. The issue we must address is whether it makes 
sense to define a perturbative S'-matrix on a time-dependent spacetime, due to the inevitable 
brmation of spacelike singularities or black holes, as predicted by the singularity theorems 



lidlnl. 



A. Colliding plane waves 

One of the simplest and best studied time-dependent geometries is the pp-wave. This is 
a solution of Einstein's vacuum equations with metric given by 

ds 2 = 2dudv + H(u, X)du 2 + dX { dX h (25) 
10 



where H(u,X) is a harmonic function V 2 H(u, X) = 0. These solutions may easily be 
extended to include nontrivial dilaton and NS 2-form sources. However, for simplicitly of 
presentation we shall concentrate on pure vacuum solutions. We will focus on the class of 
"exact" plane waves for which H(u,X) = A i j{u)X l X : > with YliAu{u) = 0. All pp-wave 
spacetimes admit a null Killing vector and this allows a global definition of null time. 
As a consequence, the natural "in" and "out" field theory vacua are chosen by decomposing 
modes into positive and negative frequency with respect to this time, and since this is a 
global deflation we find that for a free fie.d |0),„ = |0> o „« implying the absence of part.de 
creation |26j ]. (This is true provided we make the conventional choice of vacua. If the plane 
wave geometry is not asymptotically Minkowski at past and future infinity then the reasons 
for this choice are less clear). 

In order to obtain a background with pair creation we will consider a spacetime describing 
the collision of two plane waves. Treating the amplitude of each wave A^' as small, then at 
the linearized level this will be described by the metric, 



At the linearised level, the two waves pass through each other undisturbed. The 0(A 2 ) 
terms will include interactions between the waves. If the amplitude of each wave is localized 
around a given null time, for example if we take 



then at leading order the interaction region will be localized near (\u\ < L, \v \ < L). However, 
at higher orders we will see that there is a long tail in the future lightcone of the interaction 
region in which perturbations grow, eventually leading to a singularity. 

Fortunately all this can be seen explicitly in a concrete example where Af- and A^ are 
both diagonal. In this case it is straightforward to lift the linearised metric (|26|) to an 
exact nonlinear solution of Einstein's equations. To see this, first perform a linearized gauge 
transformation to take the metric to the form 



ds 2 = 2dudv + AjAu^Xidu 2 + ArJv^Xldv 2 + dX l dX t + 0{A 2 ). 



(26) 



A+Au) « Aexp(-u 2 /L 2 ) 



(27) 



ds 2 = Idudv + (S i:i + h i:j (u, v))dX i dX j + 0{A 2 ). 



(28) 



where 




(29) 



11 



This metric is a special case of a more general class of solutions given by 

ds 2 = 2e rM dudv + a{u, v) e m(u ' v) dxf , (30) 

i 

where Yli A = 0- This is the form of an exact set of solutions to Einstein's equations known 
as Einstein-Rosen waves. They satisfy d u d v a = and <9 n (a<9„A) + d u (ad v /3i) = 0. This form 
is invariant under transformations u — > f(u), v — > g(v) and exp(T) — > exp(T)/(f\ug,v)- We 
may use this freedom to set a = (u + v)\/2 = t and z = [u — v)\/2 then the satisfy 

^d t (tdA) - d% = 0. (31) 

It is straightforward to see that as t — > the solutions to this equation behave as /3j — > Ci(z) + 
di(z)\nt. The resulting solution describes a Kasner metric where the Kasner exponents 
are functions of z, consequently the generic collision of two plane waves will result in an 
anisotropic singularity. 

The existence of a singularity should cause us no surprise, it arises as a simple consequence 
of the singularity theorems. Although the singularity theorems do not strictly apply here 
as we have no matter, compactifying the geometry on one direction gives us a theory with 
gravity and a scalar field which does satisfy the singularity criteria. In the present context 
it implies that perturbation theory will inevitably breakdown in a finite time period. In 
the traditional definition of an S-matrix theory, the "out" state is defined at future infinity, 
thus it would seem at first sight impossible to define an S-matrix on a spacetime with a 
singularity at future infinity. 

B. Perturbative solution 

Our approach to dealing with this fundamental issue is a pragmatic one, for if there 
exists a long period of time for which perturbation theory is valid, then the S'-matrix should 
describe physics in that period. Let us now construct the perturbative solution to Einstein's 
equations. This will enable us to check that nonlinear corrections remain small at low orders 
in the perturbation expansion, and to explore the breakdown of perturbation theory. To be 
consistent with the metric fl2EJ), we take p i ~ 0(A), a ~ 1 + 0{A 2 ) and V ~ 0{A 2 ). The 
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perturbed metric up to and including 3rd order in perturbations is, 



a = 1 + A 2 (a + (u) + a^(v)) + 0(A*), (32) 
P = A(Pi(u) + (3i(v)) - l -A\a + {u)Pl{v) + a_(«)/3*.(u)) + 0(A 4 ), (33) 

r = -Ia 2 ^/3;(u)/^) + 0(A 4 ), (34) 



where, 



i 

oc'L{v) = -\Y.P%- (35) 

i 

We are assuming that /3+(ti) and (3-(v) are localized functions, for instance with a Gaussian 
profile /3 + = exp(— u 2 /L 2 ) and /3_ = exp(— u 2 /L 2 ). From this it is clear that up to 3rd order 
in A, /3 l and T are both bounded functions localized near the interaction region u — v — 0. 
However, as a consequence of (|33|) . a is not bounded. We choose the solution 

a+(«) = -i / M dux ^^ 2 U2 (« 2 ), (36) 

and similarly for a_(t>). From this it is clear that as u — >• +oo, then a+(w) — > C + + -D+m, 
where, 

1 f°° 

Trying to remove this divergence by adding the homogeneous solution —D + u causes a_(it) 
to diverge as u —>■ —oo. This is the first indication that perturbation theory breaks down 
at late times. On dimensional grounds, D + ~ 0(1) /L and so this coordinate system breaks 
down at u ~ L/A 2 , v ~ L/A 2 . Evidently by making A sufficiently small, we can push this 
region arbitrarily far into the causal future of the interaction region. 

Remarkably, to 3rd order in A, the growth of a is a coordinate artifact and does not 
reflect a genuine breakdown of perturbation theory. One may show this by removing the 
linear growth in u by means of the following coordinate transformation, 

A 2 

x l -> x l (l + —(D + u + D-v)) + 0(A 3 ), (38) 



A 2 

u -> u+—D + x 2 + 0(A 3 ), (39) 
A 2 

v -> v + — D„x 2 + 0(A 3 ). (40) 
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Similarly all components of the Riemann tensor are finite and supported only near u = v = 0. 
In fact we must go to 4th order in perturbations to see the first signal that perturbation 
theory is breaking down. An explicit calculation shows that certain components of the 
Riemann tensor diverge linearly for large u and v in any orthonormal frame. However, 
unsuprisingly this only occurs when u ~ L/A 2 , v ~ L/A 2 . So again we may push the 
inevitable breakdown of perturbation theory off to arbitrarily far in the future of the region 
of interest. 

In what follows we shall only compute the string amplitudes to second order in the 
background perturbations. As we have seen, at second order the collision of two plane 
waves results in a localized interaction region and no pathologies occur. Thus there will be 
no problem interpreting the implications of our amplitudes. If we continue the amplitude 
calculations to higher order in perturbations, we may expect to see mildly pathological be- 
haviour associated with the eventual breakdown in perturbation theory. However, it should 
always be possible to separate this from the information which describes the interaction 
region. These above observations show the fundamental limitations of the application of the 
S'-matrix formalism to a time- dependent geometry, these issues would not arise in a more 
Schrodinger like prescription where we consider the state at a given time, rather than the 
transition amplitude to a state at future infinity. This suggests that string field theory or a 
similar formalism may be a more appropriate way to consider time-dependent spacetimes. 



IV. PAIR PRODUCTION 



The first nonvanishing contribution to string pair production arises from the four-point 
amplitude with two incoming gravitons. In this section we discuss this amplitude for an 
explicit set of string states, and obtain results in accord with previous investigations on the 
Ugh -.„e rgy behavior of string Bering ampMudes flQQ. T his suggests that the 
behavior we observe here may hold for more general excited string states as well. A detailed 
calculation of the relevant amplitude may be found in Appendix [X] 

We calculate the pair creation of "representative" string states given in oscillator notation 

by 

|e;*> = (A!)-^,..^ I U a%at\ ] \0;k), (41) 
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where we assume that N > 1. This set of states include a unique scalar state at each 
exciation level N, which provides the most direct comparison with field theory results. 

Clearly the polarization e must be symmetric under the interchange of two holomorphic 
indices. The physical state conditions are satisfied provided that 

m 2 = -L(jV-l). (42a) 
a 

</-.,-:...,-..,- v,,v = °> for a11 3- (42b) 

'/"'"'..,,,.,„,.. ('. for all j, k. (42c) 

along with similar conditions for the barred indices. These states will have unit norm 
provided that 

(e; k\e- k) = 1 if e^^^^e^...^...^ = 1. (43) 
The vertex operators corresponding to these states are given by 

e ik - x (z,z): (44) 

where we have included a factor of the closed string coupling g c as is conventional. 

As we are interested in comparing string results to those in field theory, we focus on the 
production of "long wavelength" string states, or equivalently those whose spatial momenta 
are small in comparison with l/\fa'. Since we are also considering the production of very 
massive states, the relevant process is therefore one in which the pair of strings are created 
nearly on threshold. In order to fix notation, we take k± to be the momenta of the cre- 
ated strings, and — ki,,R the momenta of the incoming gravitons, so that the equation for 
conservation of momentum is 

k + + k^ + k L + k R = 0. (45) 

We further take the created strings to be scalar representative states of excitation levels 
N±. For simplicity of exposition we take the incoming gravitons to possess momentum in 
the t, x plane only. Working in the center of mass frame, and focusing on the case where 
the outgoing strings are created nearly on threshold, it is convenient to parameterize the 



Vjv(fc) = g c e 



jUl/Xl.../ijv/ljv \ / 



N 



N 



YldX^dXto 

L?'=i 
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momenta as 



k + = (m+, 0, 0) + (5oo + , 5k x , 5k T ), (46a) 

k_ = (m_, 0, 0) + (8uo~, -5k x , -5k T ), (46b) 

k L = (-co, -co, 0) + (-5oo L , -5k' x , 0), (46c) 

k R = (-uj, +uo, 0) + (-5oo R , +5k' x , 0), (46d) 

Momentum conservation and the mass-shell conditions imply that that we may choose Sk x 
and Skx to be the independent variables in the problem. 
For our discussion, we will find it useful to define 

a't 

a = N + -—, (47a) 
b = AL - (47b) 

where t and u are the conventional Mandelstam variables. For the problem at hand, when 
the strings are produced with zero spatial momenta, these variables are given by 

s = m 2 + + m 2 ^ + 2m + rri-, (48a) 
t = —m + m_, (48b) 
u = —m + m_. (48c) 

It will be important for our argument below that in the case where the pair is produced 
with zero spatial momentum, these variables are all 4/a' multiplied by integers. 

With this parameterization of the momenta, standard techniques in string perturbation 
theory lead to the leading order term in the string 5*-matrix, given by 

e + - - e~- - e z, -e fl -(JV + n~ 1 (iV !) _1 x 

,n N++N- N. 



\ / o — 1 ; — i 

sin (tth ] sin ( ttK\ 

■ (49) 



3=1 k=l 

T(q)r(6) 

T(a + b) 



1 2 sin(7ra) sin(7r6) 
sin[7r(a + b)] 



This amplitude is then multiplied by the metric perturbations and integrated according to 
equation 1)120 . There are two key features of this amplitude that differ substantially from 
the field theory case. 
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The first remarkable feature of this amplitude is that the production of identical string 
pairs vanishes at zero spatial momentum. This follows directly from the behavior of the 
sine functions appearing in the amplitude (|49|). When N + = iV_ = N, and when the 
spatial momentum of the created pair vanishes, then using (J47|48j) and the string mass shell 
condition we find 

o = 6 = 2iV-l, (50) 

and therefore the trigonometric factors vanish. For field theory, examination of the .S-matrix 
element arising from (J7J) reveals that the pair production amplitude is independent of k x ^ 
in the long-wavelength limit. For strings, by contrast, it appears that the pair production 
vanishes as a power law in the long-wavelength limit. Similar behavior obtains whenever 
the condition a, b G Z is satisfied, which occurs for pairs of excitation levels such as 

(JV+, AL) = (2, 5), (2, 10), (2, 17), (3, 9), . . . (51) 

in addition to the N + = AL cases corresponding to the creation of identical pairs. When N± 
are such that a ^ Z and b ^ Z, then the production of these string pairs are independent of 
spatial momentum in the long wavelength limit, similar to the field theory case. 

This behavior is a consequence of some fundamental features of string theory. Recall 
that the open string Veneziano amplitude is essentially uniquely determined by the presence 
of poles corresponding to excited string states, as well as worldsheet duality. The closed 
string amplitudes are in a sense products of open string amplitudes, with additional sine 
factors to ensure that all of poles that occur in the product are simple ones. Again, the 
closed string amplitudes are essentially uniquely determined by basic features of the theory. 
Thus different sine factors are zero when the s, t or u Mandelstam variables correspond 
to a physical string state and compensate for some of the poles arising from the gamma 
functions in the amplitude. In our case, when the string pairs are produced with zero 
spatial momentum, the s Mandelstam variable corresponds to an on-shell string state, but 
the t and u variables are the negative of on-shell values. Thus the sine functions have zeros, 
but now the gamma functions have no poles, and so the amplitude vanishes. Of course, 
this argument requires that we establish exactly some additional integers appearing in the 
arguments to the gamma functions, and this is done in Appendix [XJ The vanishing of this 
amplitude holds not only for the leading term in the 5-matrix, but for terms of all orders 
in k appearing in this amplitude. 
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The second key difference between the string and field theory pair production amplitudes 
is an exponential suppression of the string pair production rate. To explore this feature we 
consider the case where excitation level of both strings is approximately N, in which case 
a ~ b ~ 2N from (50). Applying Stirling's approximation to the gamma function factors in 
(49) we find 

r(4iV) 2 ~ N 2 \ + 8N + J 7V>>1 - (52j 
Therefore, the string amplitude falls off exponentially with m 2 , in contrast to the polyno- 
mial dependence of the field theory result. This behavior is characteristic of hard scattering 
string processes. As we have remarked earlier, typically we expect pair production rates in 
field theory to fall off exponentially with mass. However, recall that when studied with the 
S'-matrix approach employed herein, there are actually two effects at play. The first is the 
structure of the field theory amplitude connecting the incoming gravitons to the outgoing 
particles. The second is the number of gravitons, or Fourier coefficients of the metric per- 
turbations. In field theory, the former behaves polynomially with mass, while the second 
drops off exponentially. For strings on the same background (that is, with the same Fourier 
coefficients for metric perturbations) both factors fall off exponentially. Thus fewer pairs 
are produced. 

Again, this may be understood using simp 
tudes, one of the kinematic regions of interest 



e strin g p hysics. In studies of string ampli- 
14 , lslllq is the so-called "hard scattering" 



limit 

s — > oo, t/s fixed. (53) 

This is precisely the limit of interest when we look at particles of the same mass and look 
at the limit in which m — > oo. In the hard scattering limit, string amplitudes behave as 

amplitude ~ e _s/(9) (54) 

where f(8) > 0. This is precisely the behavior observed in our results. It indicates that, 
just as the scattering of strings is softer than that of point particles at high energy, so too 
the pair production of strings is reduced. 
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V. CONCLUSIONS 



In this work we have studied the pair production of excited strings in a time-dependent 
background. The specific background we consider consists of two colliding plane gravita- 
tional waves. Through studying the singularity structure of this background, as well as the 
nature of corrections needed to cancel "tadpole" diagrams, we have concluded that reliable 
results may be obtained through standard S-matrix methods at second order in the grav- 
itational wave amplitude. Our calculations revealed two essential differences between the 
field and string theory results. The pair production of identical string states is suppressed 
for certain pairs of excitation numbers of the outgoing strings; most significantly, the pro- 
duction of identical string states vanishes in the infinite wavelength limit. In addition, the 
overall production of strings is suppressed relative to point-particle results, which may be 
viewed as a consequence of the mild hard scattering behavior of string amplitudes. 

Our results are especially relevant to questions regarding the pair production of strings 
in cosmological spacetimes. In particular, these results may have implications for one mo- 
tivation for the current work, in which the pair production of strings is studied using an 
effective field theory approximation [l|, 3, 3]- The suppression of production of excited string 
states suggests that, despite the exponentially growing Hagedorn density of states, effective 
field theory provides an overestimate of the production of individual string states. 

Of course, we have established these results using techniques that are reliable only when 
spacetime is nearly Minkowski. Nevertheless, our findings suggest phenomena that may also 
be present in the strong field regime. It would be quite interesting to know if these effects 
persist in more highly curved, time-dependent backgrounds. If so, it would be an example 
of a uniquely "stringy" signature of relevance to cosmology, providing further clues to the 
role of quantum gravity in the early universe. 
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APPENDIX A: gg -► NM 



Standard techniques in string perturbation theory give the form of the four-point string 
S'-matrix amplitude as 



57T2 



M = —(: c(z + )c(z + )Vt + (z + ) :: c{z^)c{zJ)V^_ (z.) :: c(z L )c(z L )V^ L (z L ) :: V«(z R ) :> 

(Al) 

where the fact that this expression is to be integrated over the unfixed vertex operator 
position is understood. We have also suppressed the momentum conservation factor, so that 
the full ^-matrix is 



S 4 (k + , fc_, k L , k R ) = (27r) 2b 6 2b (k + + k_ + k L + k R )M 



(A2) 



In this work, we take V^ + (z + ) and V^_(zJ) to be representative states at level N + and AL, 
respectively, and Vfc L (z L ) and V£ (zr) to be the massless graviton states. We fix the three 
vertex operator positions z + = 0, = 1 and zl = oo. Next, using momentum conservation 
and the transverse property of the polarization tensors of each state, it is possible to put 
the matrix element in the form 



jV 4 



j'=i 



: ][ <).\" + 



la 



k N 

K H _|_ 

z 



fc=l 



:: JN - 



~2~ 



K 

Z 



1-Z 

k T 



X 



X 



z{\-z) 

z a'k+-k R /2^ _ z y'k--k R /2^ ^g-j 



multiplied by its complex conjugate from the antiholomorphic sector, multiplied by an overall 
factor 

and of course the appropriate polarization tensor for each state. In this expression we have 
defined, n a = k a + + k?_. This expression for the amplitude results from evaluating the ghost 
path integral and contractions between the dX and e lkX terms, as well as between two e lkX 
terms. The remaining contractions are those between two dX terms. 

Given this pattern of contractions, each possible contraction of operators results in an 
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integral of the form 



\b— 1+7712 



6-l+n 2 j2 



- 27T + mi )^(fr + m 2)T(l - a - b - m - n 2 ) 

7F r(a + b + mi + m 2 )r(l - a - m)r(l -b-n 2 ) 

For m 1 ,m 2 ,n 1 ,n2 G Z. We can simplify our discussion by focusing on the leading terms in 
this expression. Our concern in the present work is on the production of long-wavelength, 
highly excited string states. Thus, we parameterize the momenta in the problem as 



k + 
k_ 

k L 
k R 



(m+,0,0) + (5w + ,5k+,5k£), 
(m_, 0, 0) + (5lu~, 5k~, Sk?), 
( — to, -u, 0) + (-8u L , -5h%, -Sk?), 
{-u, +u, 0) + (-8u R , -8k*, -8k$), 



where we assume that 5k, 8u <C m±. We focus on the case where 

= 8% = 6k$, 
5k x = 5kt = —8kl, 



(A6a) 
(A6b) 
(A6c) 
(A6d) 



(A7a) 
(A7b) 



Now applying the mass-shell conditions and momentum conservation, we find, to = m ++ m 



and furthermore, 5k T = 5k^ = —Sk T , and 5k' x = 8k x = —8k R . We can take the two 



remaining free variables to be 5k x and <5fcr- The mass-shell conditions further imply, 

± (8k x ) 2 + (8k T ) 2 



8u 



2m± 



+ 0(kt tT ) 



(A8) 



and 5lu l = 5k' x , 5u R = 5k' x . So now we can reexpress all variables in terms of 5k x and 5kr, 
by noting, 8k' x = 5w+ + Sw . Now returning to the expression (|A5|) . we can see that 



a - 1 = {a'/2)k+ ■ k R = A 4 
b - 1 = (a'/2)h- ■ k R = N. 



1 + A + 8+r, 
1 + A + 6. R , 



with A = [(N+ - 1)(JV_ - 1)] 1/2 and 



a m 



-R 



-8u + + 



m^ 



2m_ 



a I m + 
2 V 2 



5 +R = — ( -^5uj + 



m_ 



2m A 



5uj — ujdk,. 



8ui + + uj8k x 



(A9a) 
(A9b) 



(AlOa) 
(AlOb) 
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At this point we have all of the information we need to compute each term in the amplitude. 
In order to find the dominant term, it suffices to examine each of the possible contractions 
between momentum vectors and polarization tensors. These are given by, (k R ■ e + , k R ■ e~) = 
O(\fo!rri) and (n-e + , n-e~, K-e L , k_ -e L , K-e R , k_ -e R ) = 0(\fa'5k). Thus it is seen that the 
term with the greatest number of k R contracted with polarization tensors will be dominant. 
It will be convenient to rewrite the integral (jA5|) as 

Via + mi)r(6 + m 2 )r(l — a — b — n\ — n 2 ) 

7T — -. = 

T(a + b + mi + m 2 )T(l — a — rii)T(l — b — n 2 ) 

sin[7r(a + ni)] sin[7r(6 + n 2 )] T(a + mi)T(b + m 2 ) T(a + rii)r(6 + n 2 ) ^11) 
sin[7r(a + b + ri\ + n 2 )} T(a + b + m\ + m 2 ) Via + b + ri\ + n 2 ) 

From our discussion above, the leading term with have the largest possible number of kn 
contractions, implying m x = rii = —N + and m 2 = n 2 = —N_. (This depends somewhat on 
the pattern of contractions involving the graviton polarization tensors, but this turns out to 
be unimportant in the final result) Thus we find the value of the integral to be 

sin[7r(A + 5 +R )] sin[7r(A + 6- R )] 
sin[7r(2A + 5 +R + <L fl )] 

From this it is possible to work backwards and obtain the full S-matrix expression, although 
we will not require the full amplitude. 

A key feature of this amplitude is its behavior when A 6 Z. The leading term derived 
above clearly vanishes as 8± R — >• in this case, due to the structure of the zeros of the sine 
functions. While we have only displayed the leading order term above, this behavior in fact 
holds for all of the terms contributing to this amplitude. This can be seen from the structure 
of the integral over the complex plane used to derive (jA12|) . Regardless of the pattern of 
contractions corresponding to a given term, the arguments to all of the gamma functions 
appearing in (|A12|) are positive. Since the gamma function has no zeros, and only has poles 
for negative values of the argument, the amplitude must vanish when A G Z and 5±r — > 0. 



r(A + ^ +R )r(A + ^_ R ) 

T(2A + 5 +R + 5_ R ) 



(A12) 



APPENDIX B: PAIR PRODUCTION VIA THE 5 MATRIX 

In this appendix we carry out explicit calculations of the pair production rate using both 
S'-matrix and Bogolubov coefficient techniques. This provides a check that, while the two 
methods differ, the end result is precisely the same (for a related analysis, see also (27 1). We 
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study a minimally coupled scalar field in D ^ 2 spacetime dimensions, with action 

S * = ~\f (^<R0 + m 2 2 ) d D x, (Bl) 

Our gravitational background is defined by 

9ixv = aitfr]^ (B2) 

where the "scale factor" a(t) is an arbitrary function of time. In order to better make contact 
with the examples considered in this work, we will assume that a(t) is unity, except for a 
localized "bump" of small amplitude near t — 0. Thus, there will be well-defined "in" and 
"out" regions where the geometry is Minkowski space. 

For future convenience we define a field variable ip, given by 

(j)(t, Xj) = a(t) l ~ D,2 i)(t, xj), (B3) 

We Fourier expand along the spatial coordinates, so ip(t,Xj) = if}k{t)e lkjxJ , and define 
g (t) = a 1 /™, n = -^j, and 

A(t) = {g{tf n - 1) m 2 - I (B4) 

g 

With these replacements, the action for ip becomes 

1 f ( *o 0,0 » / \ j o\ d k 



Sip — — - 



2 J \ rfc fc ™ v /r V (2tt) 
where u; 2 . = m 2 + /c 2 as usual. The equation of motion for ip k is 



i>+ [4 + A(t)]^ = 0. (B6) 



1. Bogolubov Coefficient Calculation 

Now, we will calculate precisely what the f3 Bogolubov coefficient is in the theory de- 
scribed in the previous section. Because the perturbation to the action is localized in time, 
the canonical Minkowski vacuum is a natural choice for the "in" and "out" regions, and so 

Uj {t,x)= Uj {t,x)= e - luJt+lk ^ , u* {t,x)= u* (t,x) = e +luJt - lk > xJ , (B7) 

We now consider the following (approximate) solution to the equation of motion for ip 

Mt) = (e~ tuJkt + hMt)^) (B8) 
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along with the boundary condition that fik,-k{t) — > as t — > — oo. This boundary condition 
ensures that is a pure positive-frequency mode in the "in" region. Using the definitions 
above, it is clear that f3k,~k{t) at t = +00 is the Bogolubov (3k,-k coefficient. 

Substituting our solution (jB8|) into the wave equation, and making the assumption that 
(3 < 1, we find 

/V* + ^kPk,-k = -A(t)e- 2 ^* (B9) 
This equation has the solution 

&,_*(*) = - r e ~ W r A (0 - 7T^A(2^), (BIO) 

Jt'=-oo Jt"=-oa ^k 

where in the last step we have taken t —>■ +00. Note that our derivation is self-consistent 
even in A(t) is large; we have only assumed that fait) <C 1 and thus we only need require 
that A(u)/u is small. 

To compare with the 5*-matrix calculation, we must switch to the Lorentz-invariant 
particle states, defined by 

I A;) = v/2^a+|0). (Bll) 

When this is done, one arrives at the invariant matrix element describing the particle pro- 
duction process, 

M (|0) -> |ife, -A;)) = -A(2wjfe) (B12) 

We will see in the next section how precisely the same result is obtained using conventional 
field theory techniques. 

2. S matrix Calculation 

The S'-matrix element may be calculated using the rules of standard flat-space quantum 
field theory. The calculation is much simpler than using the Bogolubov coefficient method. 
In quantum field theory, if we have a perturbation to the action 8S[i/)], then 

S= I (OlitfS'tV'piAfe) (B13) 

J vol 

where we "contract" any ip appearing in the perturbation according to the rule, 

^(a?)|ibi) = e- <fcl - !e |0> (B14) 
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In our situation, we have 5S = — A(t)-?/> 2 /2. Again, after including the factor of two from 
summing over the two possible contractions with the final state momenta, we find that, 

S = -i J e- i(kl+k2) - x A(t) d D x, (B15) 

which is precisely the Fourier transform of A(i). Imposing momentum conservation, we find 
the invariant matrix element, 

M (|0) -> \k,-k)) = -A(2uj k ) (B16) 
which is precisely the one found by other means in the previous section. 
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